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Anisotropy, thermal and quantum fluctuations and their dependence on dopant concentration 
appear to be present in all cuprate superconductors, interwoven with the microscopic mechanisms re- 
sponsible for superconductivity. Here we review anisotropy, in-plane and c-axis penetration depths, 
isotope effect and conductivity measurements to reassess the universal behavior of cuprates as re- 
vealed by the doping dependence of these phenomena and of the transition temperature. 

Establishing and understanding the phase diagram of cuprate superconductors in the temperature-dopant concen- 
tration plane is one of the major challenges in condensed matter physics. Superconductivity is derived from the 
insulating and antiferromagnetic parent compounds by partial substitution of ions or by adding or removing oxygen. 
For instance La 2 Cu04 can be doped either by alkaline earth ions or oxygen to exhibit superconductivity. The em- 
pirical phase diagram of L^-zSr^CuC^ [1-9] depicted in Fig. la shows that after passing the so called underdoped 
limit (x u « 0.047), T c reaches its maximum value T c (x m ) at x m ~ 0.16. With further increase of x, T c decreases and 
finally vanishes in the overdoped limit x ~ 0.273. This phase transition line is thought to be a generic property of 
cuprate superconductors [10] and is well described by the empirical relation 



T c (x) = T c (x m ) ( 1 - 2 ( 1 



2T C (x m ) 



(x Xu) \X Xj , 



(1) 



proposed by Presland et al. [11]. Approaching the endpoints along the x-axis, L^-xSr^CuC^ undergoes at zero 
temperature doping tuned quantum phase transitions. As their nature is concerned, resistivity measurements reveal 
a quantum superconductor to insulator (QSI) transition in the underdoped limit [12-18] and in the overdoped limit 
a quantum superconductor to normal state (QSN) transition [12]. 
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FIG. 1. (a) Variation of T c for I^-zSr^CuCU. Experimental data taken from [1-9]. The solid line is Eq.(l) with 
Tc (x m ) = 39K. (b) 7t versus x for I^-xSr^CuCU. The squares are the experimental data for 7t c [1,2,4,6,7] and the tri- 
angles for 7t=o [8,9]. The solid curve and dashed lines are Eq.(2) with 7t c ,o = 2 and 7 T=0i o = 1.63. 



Another essential experimental fact is the doping dependence of the anisotropy. In tetragonal cuprates it is defined 
as the ratio 7 = £, a b/£,c of the correlation lengths parallel (£ a b) and perpendicular (£ c ) to CuC>2 layers (afr-planes). In 
the superconducting state it can also be expressed as the ratio 7 = A c / X a b of the London penetration depths due to 
supercurrents flowing perpendicular (A c ) and parallel (A Q {, ) to the a6-planes. Approaching a non-superconductor to 
superconductor transition £ diverges, while in a superconductor to non-superconductor transition A tends to infinity. 
In both cases, however, 7 remains finite as long as the system exhibits anisotropic but genuine 3D behavior. There 
are two limiting cases: 7=1 characterizes isotropic 3D- and 7 = 00 2D-critical behavior. An instructive model 
where 7 can be varied continuously is the anisotropic 2D Ising model [19]. When the coupling in the y direction 
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goes to zero, 7 = £ x /£y becomes infinite, the model reduces to the ID case, and T c vanishes. In the Ginzburg- 
Landau description of layered superconductors the anisotropy is related to the interlayer coupling. The weaker this 
coupling is, the larger 7 is. The limit 7 = 00 is attained when the bulk superconductor corresponds to a stack of 
independent slabs of thickness d s . With respect to experimental work, a considerable amount of data is available on 
the chemical composition dependence of 7. At T c it can be inferred from resistivity (7 = £, a b/£,c — yPab/pc) and 
magnetic torque measurements, while in the superconducting state it follows from magnetic torque and penetration 
depth (7 = A c /A a b) data. In Fig. lb we displayed the doping dependence of gammar evaluated at T c (7r c ) and 
T = (7t=o)- As the dopant concentration is reduced, 7t c and 7t=o increase systematically, and tend to diverge 
in the underdoped limit. Thus the temperature range where superconductivity occurs shrinks in the underdoped 
regime with increasing anisotropy. This competition between anisotropy and superconductivity raises serious doubts 
whether 2D mechanisms and models, corresponding to the limit 77^ = 00, can explain the essential observations of 
superconductivity in the cuprates. From Fig. 1 it is also seen that (x) is well described by 

1t{x) = — , (2) 

X Xu 

where 7T.0 is the quantum critical amplitude. Having also other cuprate families in mind, it is convenient to express 
the dopant concentration in terms of T c . From Eqs.(l) and(2) we obtain the correlation between T c and jt- 

Tc{x) ( 7T Qm) A 2 , s 7T,0 ,„s 
= 1 - TT ~ 1 ) 7T {Xrn) = (3) 



T c {x rn ) \ 7 T (x) 

Provided that this empirical correlation is not merely an artefact of L^-zSr^CuO^ it gives a universal perspective 
on the interplay of anisotropy and superconductivity, among the families of cuprates, characterized by T c (x m ) and 
7t (xm)- For this reason it is essential to explore its generic validity. In practice, however, there are only a few 
additional compounds, including HgBa2Cu04+a [20], for which the dopant concentration can be varied continuously 
throughout the entire doping range. It is well established, however, that the substitution of magnetic and nonmagnetic 
impurities, depress T c of cuprate superconductors very effectively [21,22]. To compare the doping and substitution 
driven variations of the anisotropy, we depicted in Fig. 2 the plot T c (x) /T c (x m ) versus 7 t {x rn ) / 7 t {%) for a variety 
of cuprate families. The collapse of the data on the parabola, which is the empirical relation (3), reveals that this 
scaling form appears to be universal. Thus, given a family of cuprate superconductors, characterized by T c (x m ) and 
It {x m ), it gives a universal perspective on the interplay between anisotropy and superconductivity. 
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FIG. 2. Tc (x) /T c (x m ) versus 7t {x m ) / jt {x) for various cuprate families: La2-^Sr^Cu04 (•, T c (x m ) = 37K, 
7Tc W = 20) [1,2,4,6,7] , (O, Tc (x m ) = 37K, 7T =o (x m ) = 14.9) [8,9], HgBa 2 Cu0 4 + a (A, T c (x m ) = 95.6K, 7Tc (x m ) = 27) 
[20], Bi 2 Sr 2 CaCu 2 8+i (★, T c (x m ) = 84.2K, 7Tc (x m ) = 133) [23], YBa 2 Cu 3 7 _5 (♦, T c (x m ) = 92.9K, 7Tc (x m ) = 8) [24], 
YBa 2 (Cui_ a Fe a ) 3 7 - 5 (■, T c {x m ) = 92.5K, 7Tc (x m ) = 9) [25], Yi-yPr^BaaCugOT-i (▼, T c {x m ) = 91K, lTc (x m ) = 9.3) 
[26], BiSraCai.yPry^Og (-4, T c {x m ) = 85.4K, 7 t=o (x m ) = 94.3) [27] and YBa 2 (Cui_ !/ Zn y ) 3 7 - S (►, T c {x m ) = 92.5K, 
7 t=o (x m ) = 9) [28]. The solid and dashed curves are Eq.(3), marking the flow from the maximum T c to QSI and QSN 
criticality, respectively. 

Close to 2D-QSI criticality various properties are not independent but related by [13-18] 

<!>§i?2 d s v 
l6ir 3 k B y ab (0) 



2 



independently of the nature of the putative quantum critical point, fee is the Boltzmann constant, <E>o the elementary 
flux quantum, X a b (0) the zero temperature in-plane penetration depth, z is the dynamic critical exponent, d s the 
thickness of the sheets, and V the correlation length critical exponent of the 2D-QSI transitions. 5 measures the distance 
from the critical point along the x axis (see Fig. la), and i?2 is a universal number. Since T c cx d s /\ 2 ab (0) cx n° , where 
n° is the aerial superfluid density, is a characteristic 2D property, it also applies to the onset of superfluidity in 
4 Hc films adsorbed on disordered substrates, where it is well confirmed [29]. A great deal of experimental work 
has also been done in cuprates on the so called Uemura plot, revealing an empirical correlation between T c and 
d s /X 2 ab (0) [30]. Approaching 2D-QS1 criticality, the data of a given family tends to fall on a straight line, consistent 
with Eq.(4). Differences in the slope reflect the family dependent value of d s , the thickness of the sheets, becoming 
independent in the 2D limit [14-18]. The relevance of d s was also confirmed in terms of the relationship between 
the isotope effect on T c and 1/A^ fc [16,31]. Moreover, together with the scaling form (4) the empirical relation (1) 
implies 2D-QSI and 3D-QSN transitions with zV = 1, while the empirical behavior of the anisotropy (Eqs.(2) and (3)) 
require 17 = 1 at the 2D-QSI criticality. Thus, the empirical correlations point to a 2D-QSI transition with z = 1 and 
77 = 1. These estimates coincide with the theoretical prediction for a 2D disordered bosonic system with long-range 
Coulomb interactions, where z = 1 and V ~ 1 [32-34]. Here the loss of superfluidity is due to the localization of the 
pairs, which ultimately drives the transition. From the scaling relation (4) it is seen that measurements of the out of 
plane penetration depth of sufficiently underdoped systems allow to estimate the dynamic critical exponent z directly, 

in terms of T c cx (l/X 2 c (0)) z/iz+2) , which follows from Eq.(4) with 7T = A c (0) /X ab (0). In Fig.3 we displayed the 
data of Hosseini et al. [35] for heavily underdoped YBa2Cu307_5 single crystals. The solid line is 1/Aj? (T = 0) — 2.44 
10~ 4 T C 3 and uncovers the consistency with the 2D-QSI scaling relation T c cx (l/X 2 c (0)) z ^ z+2 ' with z = 1. 




T (K) 

FIG. 3. 1/Ac (T = 0) vs. T c for YBa2Cu306+:r single crystals, taken from Hosseini et al. [35]. The solid line is 
l/Xl (T = 0) = 2.44 10~ 4 T C 3 and indicates the consistency with the 2D-QSI scaling relation T c cx (l/A* (0)) z/(z+2) and 2 = 1. 

Furthermore, Hosseini et al. [35] found that at low temperature 1/Aj: (T) — 1/X 2 (0) is nearly doping independent. 
Given this behavior close to the 2D-QSI transition, scaling predicts that in leading order 1/A^ (T) — 1/A;? (0) cx T 3 
holds, in good agreement with the experimental data displayed in Fig. 4. As the in-plane penetration depth is 
concerned, there is mounting experimental evidence that 1/A^ fc (T) — 1/A„ b (0) is in the limit T — > nearly doping 
independent as well [15]. In this case 2D-QSI scaling predicts in leading order 1/A^ 6 (T) — 1/A^ b (0) cx T, in agreement 
with the experimental data of numerous cuprates [15]. 
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FIG. 4. 1/X 2 C (T) - 1/Aj? (0) vs. T for Yf^CuaOe+z single crystals, with T c values w 20.2, 19.5, 18.2, 17.8, 16.4, and 15. IK 
taken from Hosseini et al. [35]. The solid curve is 1/\ 2 C (T) - 1/A|! (0) = 3.4 10~ 7 T 3 . 

We have seen that the doping tuned flow to the 2D-QSI critical point is associated with a depression of T c and 
an enhancement of jt- It implies that when the 2D-QSI transition is approached, a no vanishing T c is inevitably 
associated with an anisotropic but 3D condensation mechanism, because 77^ is finite for T c > (see Figs. 1 and 2). 
This represents a serious problem for 2D models [36] as candidates to explain superconductivity in the cuprates, and 
serves as a constraint on future work toward a complete understanding. Note that the vast majority of theoretical 
models focus on a single Cu-0 plane, i.e., on the limit of zero intracell and intercell c-axis coupling. 

Since Eq.(4) is universal, it also implies that the changes AT C , Ad s and A (l/A^ b (T = 0)), induced by pressure or 
isotope exchange are not independent, but related by 

AT C _ Ad A (1/X 2 ab (0)) _ Ad s A(X ab (0)) 

T c d s + (1/A^(0)) d s A a6 (0) ■ [ ' 

In particular, for the oxygen isotope effect ( 16 vs. 18 0) of a physical quantity X the relative isotope shift is defined 
as AX/X = ( 18 X - 16 X)/ 18 X. In Fig. 5 we show the data for the oxygen isotope effect in La 2 - a; Sr x Cu04 [37,38], 
Yi_ a Pr x Ba2Cu307_5 [38-40] and YBa 2 Cu307_,5 [38,41], extending from the underdoped to the optimally doped 
regime, in terms of A (A afc (0)) / A ab (0) vs. AT C /T C . 
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FIG. 5. Data for the oxygen isotope effect in underdoped La 2 - a! Sr a: Cu04(0 : x=0.15 [38], A:x=0.08, 0.086 [37], 
Yi^Pr^BaaCuaOT-s (•: x=0, 0.2, 0.3, 0.4) [38-40] and YBa 2 Cu 3 7 -i (T [38], ■ [41]) in terms of A (X ab (0)) / X ab (0) 
vs. -ATc/Tc. The solid line indicates the flow to 2D-QSI criticality and provides with Eq.(5) an estimate for the oxygen isotope 
effect on d s , namely Ad s /d s = 3.3(4)%. 

It is evident that there is a correlation between the isotope effect on T c and X a b (0) which appears to be universal 
for all cuprate families. Indeed, the solid line indicates the flow to the 2D-QSI transition and provides with Eq.(5) 
an estimate for the oxygen isotope effect on d s , namely Ad s /d s — 3.3(4)%. Approaching optimum doping, this 
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contribution renders the isotope effect on T c considerably smaller than that on A a &(0). Indeed, even in nearly 
optimally doped YBa2Cu30y_5, where AT C /T C = —0.26(5)%, a substantial isotope effect on the in-plane penetration 
depth, AA a h (0) /\ a b (0) = —2.8(1.0)%, has been established by direct observation, using the novel low-energy muon- 
spin rotation technique [41]. Note that these findings have been obtained using various experimental techniques on 
powders, thin films and single crystals. 

In this context it is important to recognize that the substantial isotope effect on the in-plane penetration depth 
is incompatible with the Migdal-Eliashbcrg (ME) [42] of the electron-phonon interaction [43]. Since the lattice 
parameters remain essentially unaffected [44,45] by isotope exchange, while the dynamics associated with the mass 
of the respective ions is modified, the substantial isotope effect on the zero temperature penetration depth requires 
a renormalization of the normal state Fermi velocity vjr — ► vp = vp/(l + /) where vp is the bare velocity and 
/ the electron-phonon coupling constant. However, in the Migdal-Eliashberg [42] treatment of the electron-phonon 
interaction / is independent of the ionic masses and assumed to be small [46,47]. This is true if the parameter 
u>of / Ep is small, where loq is the relevant phonon frequency and Ep the Fermi energy. Thus the isotope effect on the 
penetration depth is expected to be small, of the order of the adiabatic parameter 7 = uoo/Ep « 1. The ME theory 
retains terms only of order 0. Cuprates, however, have Fermi energies much smaller than those of conventional metals 
[48] so that 7 is no longer negligible small. In any case the substantial isotope effect on the in-plane penetration depth 
uncovers the coupling between local lattice distortions and superfluidity and the failure of the Migdal-Eliashberg (ME) 
treatment of the electron-phonon interaction, predicting 1/A 2 (0) to be independent of the ionic masses [42]. Evidence 
for this coupling emerges from the oxygen isotope effect on d s , the thickness of the superconducting sheets, while 
the lattice parameters remain unaffected. Indeed, the relative shift, Ad s /d s w 3.3(4)%, apparent in Fig. 5, implies 
local oxygen distortions which do not modify the lattice parameters but are coupled to the superfiuid. Although 
the majority opinion on the mechanism of superconductivity in the cuprates is that it occurs via a purely electronic 
mechanism involving spin excitations, and lattice degrees of freedom are supposed to be irrelevant, the relative isotope 
shift Ad s /d s rts 3.3(4)% uncovers clearly the existence and relevance of the coupling between the superfiuid and local 
lattice distortions. 

Finally we turn to the finite temperature critical behavior. Close to the critical temperature T c of the superconductor 
to normal state transition, in a regime roughly given by the Ginzburg criterion [14,49-51], order parameter fluctuations 
dominate critical properties. In recent years, the effect of the charge of the superconducting order parameter in this 
regime in three dimensions has been studied extensively [52-59,61]. While for strong type-I materials, the coupling of 
the order parameter to transverse gauge field fluctuations is expected to render the transition first order [53] , it is well- 
established that strong type-II materials should exhibit a continuous phase transition, and that sufficiently close to 
T c , the charge of the order parameter is relevant [55-61]. However, in cuprate superconductors within the fluctuation 
dominated regime, the region close to T c , where the system crosses over to the regime of charged fluctuations turns 
out to be too narrow to access. For instance, optimally doped YBa2Cu307_a, while possessing an extended regime 
of critical fluctuations, is too strongly type-II to observe charged critical fluctuations [14,49-51,62]. Indeed, the 
effective dimcnsionless charge e = £/A = 1/k is small in strongly type II superconductors (k >> 1). The crossover 
upon approaching T c is thus initially to the critical regime of a weakly charged superfiuid where the fluctuations of 
the order parameter are essentially those of an uncharged superfiuid or XY- model [49]. Furthermore, there is the 
inhomogencity induced finite size effect which renders the asymptotic critical regime unattainable [63-65]. However, 
undcrdoped cuprates could open a window onto this new regime because n a b is expected to become rather small. 
As outlined above, in this regime cuprates undergo a quantum superconductor to insulator transition [14-18] and 
correspond to a 2D disordered bosonic system with long-range coulomb interactions. Close to this quantum transition 
T c , A a b and £ a & scale as T c oc A~ b 2 oc £~> 2 oc A c 2 ( z+2 )/ z [14-18], while £ c — > d s . These relations yield with the dynamic 
critical exponent z = 1 [14,18,32-34], n a b oc T^ 2 and n c oc T c -3 ^ 2 . Noting that T c decreases by approaching the 
underdoped limit, the in-plane penetration depth appears to be a potential candidate to observe charged criticality in 
sufficiently homogeneous and underdoped cuprates, while the c-axis counterpart is expected to exhibit neutral critical 
behavior. 

When charged fluctuations dominate the in-plane penetration depth and the correlation length are related by [57-61] 

A a 6 = ^abiab, Kb = Kab \t\ " , V ~ 2/3, (6) 

contrary to the uncharged case, where A oc and with that 

Kb = A 0a6 \tr /2 , (7) 

where t = T/T c — 1. In a plot (din Kb/dT)^ 1 vs. T charged critical behavior is then uncovered in terms of a 
temperature range where the data falls on a line with slope \ jv ~ 3/2, while in the neutral (3D-XY) case it collapses 
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on a line with slope 2/v ~ 3. Clearly, in an inhomogeneous system the phase transition is rounded and (d In A a b/dT) _1 
does not vanish at T c . In Fig. 6 we displayed (T c dln Xab/dT)^ 1 vs. t for YBa2Cu30g.59 with applied hydrostatic 
pressures of 0.5kbar and 10.5kbar, derived from the measured in-plane penetration depth data of ref. [66]. The 
solid line with slope 1/v ~ 3/2 indicates according to Eq. (6) the charged criticality in homogeneous systems with 
T c = 57.82K and T c — 61. IK. Although the charged criticality is attained there is no sharp transition because X a b 
does not diverge at T c . Inhomogeneities prevent the correlation length £ ;, = X a b/^ab to grow beyond the spatial 
extent L a b of the homogenous domains in the a6-plane. For a discussion of the associated finite size effect we refer to 
ref. [66] . In any case these measurements clearly reveal that the temperature dependence of the in- plane penetration 
depth opens a window onto charged criticality at finite temperature in underdoped cuprates. 
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FIG. 6. (TcdlnXat/dT)- 1 with X ab in /im vs. t for YBa2Cu3 06.59 under hydrostatic pressure of 0.5kbar (A)and 10.5kbar 
(O) taken from ref. [66]. The straight line with slope 1/v ~ 3/2 corresponds according to Eq. (6) to charged critical behavior 
with T c = 57.82K and T c = 61. IK. 



Extended c-axis penetration depth measurements on Yl^CuaOg+z single crystals with T c 's ranging from 4 to 20K 
have been performed by Hosseini et al. [35]. In Fig. 6 we displayed their data for the sample with T c w 19.75K in 
terms A~ 2 and (din X~ 2 /dT)~ 1 vs. T. The solid curve is A" 2 = A 0c 2 (1 - T/Tcf with T c = 19.75K and A^ 2 = 2.15 
(10 9 m- 2 ) and the dashed one (din X^/dT)- 1 = -l/v(T c - T) with v ~ 2/3. 




T(K) 

FIG. 7. A^ 2 and (din X~ 2 /dT)~ 1 vs. T for a YBaaCusOe+x single crystal derived from the data of Hosseini et al. 
[35]. The solid curve is X~ 2 = A ( 7 C 2 (1 - T/Tc)" with T c = 19.75K and Ag c 2 = 2.15 (10 9 m~ 2 ) and the dashed one 
(dlnA^ 2 /dT) _1 = -l/v(T c — T) with v ~ 2/3. They indicate according to Eq. (7) the leading uncharged (3D-XY) criti- 
cal behavior. 

They indicate according to Eq. (7) the leading uncharged (3D-XY) critical behavior. Noting that in this underdoped 
regime the anisotropy 7 = X c /X a b is rather large (see Fig. 2) the critical regime where thermal 3D fluctuations 
dominate will be compared to optimum doping much narrower. Here it extends to t 77/93.7 — 1 ~ —0.18 [62]. 
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Nevertheless, a glance to Fig. 7 shows that the 3D-XY critical regime is attained for 18K< T < 19K, corresponding to 
0.09 > \t\ > 0.04. The tail above 19K is attributable to an inhomogeneity induced finite size effect. Indeed, according 
to the scaling relation k c a T c ' it is unlikely that the crossover to charged criticality is attainable. In contrast and 

1/2 

consistent with n a b oc T c charged criticality turned out to be accessible in the a6-plane penetration depth. 

We have seen that the linear relationship between T c and 1/Aj? 6 (0) in the underdoped regime, referred to as the 
Uemura relation, is a consequence of the the dimensional crossover associated with the flow to the 2D-QSI transition 
(see Eq.(4). To illustrate this behavior we displayed in Fig. 8 T c vs. 1/A^ fc (0) for I^-^Sr^CuO-i. The straight line 
is Eq.(4) in terms of T c = 25A^ b 2 (0) and the arrow indicates the flow to the 2D-QSI criticality. However, approaching 
the optimally doped regime, where T c reaches its maximum value (see Fig. 1), the observed T c 's are systematically 
lower than the 2D-QSI line and this trend continues in the overdoped regime. 




x ab " 2 (0) 10 7 (A" 2 ) 

FIG. 8. T c vs. l/\ab (0) for L^-xSr^CuCU. Data taken from Uemura et al. [30,67] and Panagopoulos et al. [9]. The straight 
line is Eq.(4) with R.2d s = 6.5A and the arrow indicates the flow to 2D-QSI transition criticality. 



To provide an understanding we invoke 3D-XY universality expected to hold along the phase transition line, T c (x) 
in La2- x Sr x Cu04 (see Fig. I), as long as the charge of the pairs can be neglected. In this case the transition 
temperature T c and the critical amplitudes of the penetration depths A a f,o and transverse correlation lengths £*£o are 
related by [14,18] 



167r KbO 167r A c0 



where A? (T) = X% (I - T/T c ) v and $ (T) = ^ (I - T/T c ) v with v ~ 2/3. For our purpose it is convenient to 
express the transverse correlation length to the correlation lengths above T c in terms of [14,18,68] 



where the anisotropy is given by 



^ = / » 0.453, p- - 7 /, (9) 

?c0 t,abO 



2 _ I AeO A _ CcO _ ( CabO 
v AahO / CabO V £< 



7^=(t^) = ^ = ( ™ ) • (10) 



Combining Eqs.(8) and (9) we obtain the universal relation 

_3 

I67T 3 fcs 



$ 2 f 

T c X lbO = 1c Jji. (11) 



It holds, as long as cuprates fall into the 3D-XY universality class, irrespective of the doping dependence of T c , A^ h0 
and £ c0 . Fo this reason it provides a sound basis for universal plots. However, there is the serious drawback that 
reliable experimental estimates for T c and the critical amplitudes A a fco and £ c o measured on the same sample are not 
yet available. Nevertheless, some progress can be made by noting that in the underdoped regime, approaching the 
2D-QSI transition the universal scaling forms (4) and (II) should match. This requires 
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so that away from 2D-QSI criticality 



/£ c0 -> iZd a , (12) 



«J0) -g-£ c0 , (13) 



167T 3 fcs 



holds. Thus, when both T c and 1/A^, (0) increase, T c values below T c = ( y <&lR 2 /6ir 3 k B ) d s /X 2 ab (0) (Eq.(4)) require 
£ c o to fall off from its limiting value £ c0 = {R2/ f)d s . The doping dependence of £ c o in La2_ a; Sr ;E Cu04, deduced 
from Eq.(13) and the experimental data for T c and \ 2 ab (0) is displayed in Fig. 9a in terms of T c vs. £ c o and £ c o vs. 
x . Approaching the underdoped limit (x ~ 0.05), where T c vanishes (Fig.l) and the 2D-QSI transition occurs, £ c o 
increases nearly linearly with decreasing x to approach a fixed value. Indeed, the data is consistent with 

£ c0 = (I6^ 3 fc s / ($§/)) T c X 2 ab (0) « 14.34 - 60.47(x - 0.05)A, (14) 

yielding the limiting value £ c0 (x — 0.05) w 14.34A and /£ c0 (a; = 0.05) = i?2^ s ~ 6.5A used in Fig. 8. An essential 
result is that £ c o adopts in the underdoped limit (x ~ 0.05) where the 2D-QSI transition occurs and T c vanishes 
its maximum value £ c o ~ 14.34A, which is close to the c-axis lattice constant c ~ 13.29A. Thus, a finite transition 
temperature requires a reduction of £ c o, well described over an unexpectedly large doping range by Eq.(14). Noting 
that this relation transforms with Eq.(2) to 

Ceo w 14.34 - 6O.4770/7A, (15) 

this behavior is intimately connected to the doping dependence of the anisotropy. Hence, a finite T c requires un- 
avoidably a finite anisotropy 7. The lesson is, in agreement with the empirical relation (3), that superconductivity in 
La2_a;Sr x Cu04 is an anisotropic but 3D phenomenon which disappears in the 2D limit. From the plot 1/A^ 6 (0) vs. 
T c /£ c displayed in Fig. 9b, where according to Eq.(13) universal behavior is expected to occur, the data is seen to fall 
rather well on a straight line. This is significant, as moderately underdoped, optimally and overdoped L^-xSr^CuC^ 
falling according to Fig. 8 well off the 2D-QSI behavior T c cx 1/A;? 6 (0) now scale nearly onto a single line. Thus the 
approximate 3D-XY scaling relation (13), together with the empirical doping dependence of the c-axis correlation 
length (Eqs.(14) and (15)) are consistent with the available experimental data for La2_ x Sr x Cu04 and uncovers the 
relevance of the anisotropy. However, the linear doping dependence of £ c is not expected to hold closer to the over- 
doped limit (x w 0.27) where T c vanishes and a 3D quantum superconductor to normal state (3D-QSN) transition is 
expected to occur (see Fig.l). 




Since Eq.(ll) is universal, it should hold for all cuprates falling in the accessible critical regime into the 3D-XY 
universality class, irrespective of the doping dependence of T c , X 2 ab0 , iabo, 7 and £ c o. Having seen that charged 
criticality is accessible in the heavily underdoped regime only, Eq.ll) rewritten in the form 



><■() 



£a&07 167T 3 fc S 



8 



provides a sound basis for universal plots. However, as aforementioned there is the serious drawback that reliable 
experimental estimates for the the critical amplitudes and the anisotropy at T c measured on the same sample and for 
a variety of cuprates are not yet available. Nevertheless, as in the case of L^-xSr^CuC^ outlined above, progress can 
be made by invoking the approximate scaling form (13) and by expressing the correlation lengths in terms of erf c the 
real part of the frequency dependent conductivity erf c (ui) in direction i extrapolated to zero frequency at T > T c [69] 
in terms of 

1/& * otthab, 1/ {lU) - O^/Sc, (17) 

Si with dimension reincorporates the temperature dependence. With that the universal relation (16) transforms 
with Eq.(13) to 

1 16A (18) 



K b ($)Tc<J d a i A 2 (0)T c af & fs ab 

because erf^/erf — j 2 and with that s a b — s c . To check this relation we note that the 2D-QSI scaling from (4) 
transforms with the relation d s = cr sheet /afl between sheet conductivity a sheet and conductivity afl to 

1 _ 167T 3 /c B .fceet_ h , « m-4n-l 



= —ra ~ct 1.55 HTHr 1 , (19) 



where h/Ae 2 = 6.45kSl is the quantum of resistance and ero is a dimensionless constant of order unity [72]. Thus, 

with A a b (0) in /im, T in K and cr^jj in (Q _1 cm _1 ), and the structure of the a6-expression in Eq.(18) is recovered. 
Similarly, approaching 2D-QSI criticality A c (0) and af c scale as [70] 

A c (0Hf7 s (af)- (2+z)/4 , (21) 

where z = 1 is the dynamic critical exponent of the 2D-QSI transition. Q s is a non-universal factor of proportionality. 

Noting that in this limit erf c and T c scale as T c cx (cr^ c c-expression in Eq.(18) is readily recovered. Accordingly, 

there is no contradiction between the scaling forms (18) and (21). In this context it should be kept in mind that 

the universal relation is valid for the neutral case only. However we have seen that the effective dimensionless charge 
1/2 

£ab = 1/ K a b scales as e a b oc T c and the charged critical regime becomes accessible. This is not the case for the c-axis 
penetration depth because e c oc T c 3 ^ 2 . For this reason, as the underdoped limit (2D-QSI transition) is approached 
there should be a crossover from Eq.(18) to the universal relation (19), becoming manifest in different constants on 
the right hand side. 

On the other hand, approaching the 3D quantum superconductor to normal state (QSN) transition \ a b,c (0), T c and 
erfl c scale as [18] 1/A^ bc (0) cx Tc 1+Z ^ z , afl cx y^ 2 ^ -1 )/ 2 ^ w here z is the dynamic critical exponent of this quantum 
transition and z c i the exponent associated with the finite temperature critical dynamics. Indeed, in this limit the 
correlation length £ r associated with the finite temperature critical dynamics cannot be eliminated. Since £ T scales 
as £r oc Cb we obtain af cx UrKlb « ilU) « Ct' 1 « T - {z ^ 1)/z . Accordingly, 

X 2 ab (0) T c a d a t cx (0) T c af cx (22) 

should hold close to 3D-QSN criticality. Furthermore in this regime the effective charge becomes negligible small 
so that 3D-XY scaling applies. Indeed e a b,c — ^/nab.c scales as e a b,c oc j^ 2-1 )/ 22 because X a b,c oc j , ( 7( z+1 )/ 2z an( j 
Ca6,c oc T c . A potential candidate for the 3D-QSN transition is the model proposed by Herbut [73]. It describes 
the disordered d-wave superconductor to disordered metal transition at weak coupling and is characterized by the 
dynamic critical exponents z = 2. When this holds true the effective charge scales at e a b,c oc T c 3//4 and 3D-XY scaling 
is no longer applicable. Furthermore there is experimental evidence for z c i = 2 [51,74]. Accordingly, the scaling form 
(18) does not hold in the overdoped regime close to 3D-QSN criticality. In this regime Eq.(21) is replaced by [70] 

A c (0)oc (af)^T, (23) 
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with a non-universal factor of proportionality. What is particularly remarkable is then that the 3D-XY scaling from 
(18) predicts that all points of l/X 2 ab (0) vs. T c a d a % and (0) vs. T c a dc should fall onto a single line with the 
exception of the overdoped regime where the scaling form (22) applies. 

In Fig.10 we displayed 1/A| 6 (0) vs. T c a d a l and 1/A^(0) vs. T c af as collected by by Homes et al. [71]. In 
agreement with the scaling from (18) the afr-plane and c-axis data appears to be well described by the same line, 
namely l/\\ c (0) ~ 5.2 10" 5 T c cr^ c , yielding for R 2 cr in the 2D-QSI relation (Eq.(20)) the estimate R 2 cr = 1.98. 
Furthermore, all points of 1/A^ b (0) vs. T c a d £ (open symbols) nearly fall onto a single line with a slope of unity. 
This is significant, as moderately underdoped, optimally and overdoped materials, which fell well off of the 2D-QSI 
behavior T c oc 1/A^ b (0) (see Eq.(4) and Fig. 8) now scale nearly onto a single line, in agreement with Fig. 9b. 



FIG. 10. 1/X 2 ab (0) vs. T C (T% (open symbols) and (0) vs. T c a* c (full symbols) for various cuprates as collected by Homes 
et al. [71]. The straight line is l/\% b>c (0) = 5.2 10" 5 T c cr^ c . The experimental data is taken from [75-80] for YBa 2 Cu 3 7 -i 
(Y-123), [77] for Pr/Pb-YBa 2 Cu 3 7 -,5(Pr/Pb-123), [78,81] for YBa 2 Cu 4 08 (Y-124), [78,81] for Tl 2 Ba 2 Cu0 6+ ,5, [77,82] for 
Bi 2 Ca 2 SrCu 2 8+ a and Y/Pb Bi 2 Ca 2 SrCu 2 8+i , [83,84] for Ndi.gsCeo.isCuCU, [71,85] for La^Sr^CuO^M), and [71] for 
HgBa 2 Cu04+^. 

This agreement demonstrates that the approximate scaling relation 18 captures the essentials of the exact 3D-XY 
scaling form (16) by eliminating the doping dependence of the critical amplitudes towards the 2D-QSI transition in 
terms of there zero temperature counterparts. This evidence for anisotropic 3D-XY scaling raises again serious doubts 
that 2D models [36] are potential candidates to explain superconductivity in the cuprates. However the data does not 
extend sufficiently close to 2D-QSI- and 3D-QSN criticality in order to uncover the aforementioned crossovers from 
the scaling from (18) to (19) and from (18) to (22). 

In contrast, the plot A c (0) vs. a dc displayed in Fig. 11 provides according to the scaling form (21) information 
on the flow to 2D-QSI criticality. The straight line is Eq.(21) with f) s = 24 appropriate for La2- a; Sr a ;Cu04(214) and 
the arrow indicate the direction of this flow. The data for La2-a;Sr 2; Cu04 covers the range from x = 0.08 to 0.2. 
Although the data is sparse one anticipates that the deviations from the straight line behavior increase systematically 
as the overdoped region is approached. As noted previously [70] this behavior is attributable to the initial crossover 
to 3D-QSN-criticality, where the scaling form (23) applies. Clearly more experimental data extending much closer to 
the overdoped limit are needed to uncover this crossover. Actually this also applies to the crossover from the scaling 
form (18) to (22). 




10 % 
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FIG. 11. 1/A C (0) vs. o d c c using the data shown in Fig. 10 and taken from [86,87] for YBa 2 Cu 3 C>7-i (Y-123), 
La2- a; Sr 2 :Cu04(214), and from [88] for Bi2Ca2SrCu20s+«. The straight line corresponds to Eq.(21) with Q s — 24 appro- 
priate for La2- a; Sr a; Cu04(214). 

Since Eq.(18), with the exception of the heavily overdoped regime, turned out to be nearly universal, it also implies 
that the changes AT C , A fl/A 2 bc J anc ^ A°"o6c> induced by pressure or isotope exchange are not independent, but 
related by 



^(lA^ e (0) 
(l/A 2 ofe , c (0)) 



AlT ab,c 

dc 
a ab,c 



(24) 



In particular, for the oxygen isotope effect ( 16 vs. 18 0) of a physical quantity X the relative isotope shift is defined 
as AX/X = ( 18 X - w X)/ 1S X. Since close to 2D-QSI criticality d s = a sheet /(j% and with that Aa^/a^ = -Ad s /d s 
we recover the universal relation (5). Close to optimum doping where in YBa2Cu3G , 7_ ( 5 AT C /T C = —0.26(5)% 

and A (l/A^.c (0)) / (l/A 2 6jC (0)\ = -5.6(2.0)% [41] it predicts -5.6% effect on Acrfg/cxfg evaluated above and 
extrapolated to T c . Hence, the absence of a substantial isotope effect on the transition temperature of cuprates does 
not rule out phonons as the bosons responsible for the coupling between the charge carriers, as previously suggexted 
[89]- 

To summarize we observed remarkable agreement between the experimental data of underdoped cuprates and the 
flow to the 2D-QSI critical endpoint in the underdoped limit. Although an inhomogeneity induced finite size effect 
makes the asymptotic critical regime unattainable, there is considerable evidence that the 2D-QSI transition falls 
into the universality class of a 2D disordered bosonic system with long-range Coulomb interactions. The loss of 
superfluidity is due to the localization of the pairs, which ultimately drives the transition. Important implications 
include: (i) A finite transition temperature and superfluid density in the ground state are unalterably linked to a finite 
anisotropy. This finding raises serious doubts that 2D models [36] are potential candidates to explain superconductivity 
in cuprates. (ii) The doping dependence of the large oxygen isotope effects on the zero temperature in-plane penetration 
depth confirms the flow to 2D-QSI criticality, poses a fundamental challenge to this understanding and calls for a 
theory that goes beyond Migdal-Eliashberg. Although the majority opinion on the mechanism of superconductivity 
in the cuprates is that it occurs via a purely electronic mechanism involving spin excitations, and lattice degrees 
of freedom are supposed to be irrelevant, the relative isotope shift of the thickness of the superconducting sheets 
Ad s /d s 3.3(4)% uncovers clearly the existence and relevance of the coupling between the superfluid and local 
lattice distortions. (iii) Given the striking feature that at low temperature the absolute change in 1/X ab and 1/A 2 has 
essentially no doping dependence, in contrast to l/A 2 ^ (0) and 1/A 2 , (0), scaling around the 2D-QSI critical endpoint 
implies a T 3 power law for 1/A 2 , (T) and T for l/A 2 ^ (T), consistent with experiment. Such power law behavior also 
follows from the nodes characteristic of a d-wave energy gap in the one particle density of states. However, our 
derivation shows that in underdoped cuprates these power laws are characteristic features of the 2D-QSI transition 
at the critical endpoint and hold for both s- and d-wave pairing . (iv) Similarly, we have seen that the large changes 
of l/\ 2 ab (0), 1/X 2 C (0) and T c , consistent with the relation T c cx A~ h 2 (T = 0) cx A~ fc 2/3 (T = 0), follow from the flow to 



1/2 

the 2D-QSI criticality as well, (v) Because the Ginzburg parameters scale Tc<* and K 

was shown to open a door to observe charged criticality in the in-plane penetration depth at finite temperature. An 



oc T c 3 ^ 2 this flow 
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intriguing implications for microscopic models is then the increasing value of the effective charge of the pairs in the ab- 
plane, where e a b = l/^ab ^ T c 1//2 , while in the c-direction it disappears as e c = 1/k c oc T c 3 ^ 2 . vi) We demonstrated 
here that the empirical relation \ 2 ab (0) T c a%% — \ 2 C (0) T c af c ~ const [71] is a consequence of 3D-XY universality 
extended to anisotropic systems. Nevertheless, these relations clearly reveal that the absence of a substantial isotope 
effect on the transition temperature of cuprates does not ruled out phonons as the bosons responsible for coupling 
between the charge carriers, as previously suggested. 
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